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We investigate stability of the D-dimensional Reissner-Nordstrom-anti-de-Sitter metrics as solu- 
tions of the Einstein-Maxwell equations. We have shown that asymptotically anti-de Sitter black 
holes are dynamically stable for all values of charge and anti-de Sitter radius in D = 5, 6 ... 11 dimen- 
sional space-times. This does not contradict to dynamical instability of RN-AdS black holes found 
by Gubser in A/" = 8 gauged supergravity, because the latter instability comes from the tachyon 
mode of the scalar field, coupled to the system. Asymptotically AdS black holes are known to be 
thermodynamically unstable for some region of parameters, yet, as we have shown here, they are 
stable against gravitational perturbations. 



PACS numbers: 04.30.Nk,04.50.-|-h 



I. INTRODUCTION 

Nowadays, brane world theories and string theory, im- 
ply existence of extra dimensions in nature This in- 
duced great interest to higher dimensional black holes, in 
particular, to such features as perturbations, dynamical 
and thermodynamical stability, particles and fields be- 
havior and Hawking radiation around these black holes. 
Unlike four dimensional case, higher dimensional space- 
times, admit plenty of "black" solutions: black holes, 
black strings and branes, black rings and Saturns. Sta- 
bility of these solutions may be criteria of their exis- 
tence. The stability analysis usually requires to trans- 
form the perturbed Einstein equation to the wave-like 
form, what was done yet in 1957 by Regge and Wheeler 
for = 4 black holes, and only in 2003 for general num- 
ber of space-time dimensions [2]. Then, the stability 
of D-dimensional Schwarzschild black holes was proved 
in 0] and of Schwarzschild-de Sitter black holes in 
Recently the stability of Kaluza-Klein black holes with 
squashed horizons was shown in Unlike Kaluza- 

Klein black holes, black strings and branes become un- 
stable for perturbations with wavelength, which is larger 
than some threshold value (Gregory-Laflamme instabil- 
ity [3]) ■ In [1] it was shown that the instability thresh- 
old point corresponds to some dominating static solu- 
tion of the wave equation. The neutral D-dimensional 
black holes in Gauss-Bonnet theory are unstable only for 
D = 5, 6 and for small values of Gauss-Bonnet coupling 
0j [13] ■ The instability in the Gauss-Bonnet theory is 
qualitatively different from a black string instability: the 
black string instability is an example of instability, de- 
veloped at the lowest multipoles, therefore with a static 
solution dominance at the threshold point. The insta- 
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bility in the Gauss-Bonnet theory is developed at large 
multipoles, so that the growing mode dominates after a 
long period of damped quasinormal oscillations [To| . 

The stability of the higher-dimensional black holes is 
important also for the growing interest to the quasinor- 
mal modes of the Standard Model fields in higher di- 
mensional theories [ll|. Indeed, only stable black holes 
can be considered as a background on which, test fields 
propagate. 

Higher dimensional black holes in asymptotically anti- 
de Sitter (AdS) space-times have been in the focus of 
string theorists recent years, because of their role in 
the AdS/CFT correspondence. A large asymptotically 
anti-de Sitter black hole corresponds to a thermal state 
in the dual conformal field theory, where the Hawking 
temperature of the black hole is the temperature in the 
dual field theory [l^. The perturbations of AdS black 
holes have been extensively studied during recent decade 
[T^ ]. Nevertheless, it was not known until the present 
study, if D-dimensional asymptotically AdS black holes 
are dynamically stable as solutions of D-dimensional 
Einstein-Maxwell equations. The stability of £> = 4 and 
D = 5 Reissner-Nordstrom-AdS black holes were stud- 
ied by Gubser and Mitra in the TV = 8 gauged super- 
gravity theory [l3| i -e. a theory with the Maxwell 
and scalar matter fields coupled to the electromagnetic 
field. There it was shown that the highly charged black 
holes are unstable and the parameter region of instabil- 
ity increases for larger black holes, i.e. for black holes, 
which radius is much larger than the anti-de Sitter ra- 
dius. Yet that instability evidently came from the tachy- 
onic mode of the scalar field. Thus the question remains 
if there is an instability of Reissner-Nordstrom-anti-de 
Sitter black holes within the ordinary Einstein-Maxwell 
theory? Our main aim here is to answer this question, 
keeping in mind such an important feature of black holes 
as thermodynamic (in)stability. Thus, accor ding to hy- 
pothesis of Gubser and Mitra in and [1^1, there 
may be a correlation between thermodynamic and dy- 
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namic (gravitational) (in)stabilities of black holes, be- 
cause it was found, that the parametric region of thermo- 
dynamic and dynamic (in)stabilities although do not co- 
incide, differ from each other only slightly for the Af — 8 
gauged super-gravity. Another possible correlation could 
give the thermodynamic instability of small AdS black 
holes, what happens with a phase transition, called the 
Hawking-Page transition (l6| . One of the thermodynam- 
ically preferred final states, after the transition, might be 
a pure AdS space-time [la] ■ 

The thermodynamic instability of black holes takes 
place also for ordinary Einstein-Maxwell- AdS black holes 
UM for some values of black hole parameters. The sec- 
ond order phase transition occurs at the instability point 
If we expect some correlations between thermo- 
dynamic and dynamic instabilities to these cases, we 
should expect gravitational instability of D-dimensional 
charged black holes in AdS space-times within the stan- 
dard Einstein-Maxwell theory. In this paper, we shall 
show that this is not the case of the pure RNAdS black 
holes, which are stable against gravitational perturba- 
tions for D = 5, 6 ... 11, where D is the number of space- 
time dimensions. 

The paper is organized as follows: Section II intro- 
duces the basic formula for the background metric and 
for perturbation equations reduced to a wave-like form. 
Section III reviews the numerical method, which we used 
for stability analysis. Sections IV and V consider the 
obtained results for Reissner-Nordstrom -anti- de Sitter 
black holes. 



II. BASIC FORMULAE 

The metric of the D — d + 2-dimensional Reissner- 
Nordstrom- (anti)-de-Sitter black holes is given by the line 
element 



ds^ = f{T)df - ^ - r^dn,, 



(1) 



where dQd is the line element on a unit d-sphere, and 

f{r)^l^X + Z-Y, (2) 



d{d+l) 



2d-2 ■ 



Here M is the mass parameter of the black hole, Q is 
its charge, and the A-term coincides with a cosmologi- 
cal constant, when positive, and is related to the anti-de 
Sitter radius, when negative, in the following way. 
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1. 



The general perturbations of the Einstein-Maxwell 
equations 




FIG. 1: Effective potential for scalar type of gravitational 
perturbations of the RNAdS black holes for D — 5, r+ = lOR. 
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Tg^,J) + -j^^Sg^,^, (4) 



after using the gauge freedom, and separating the angular 
variables can be reduced to a number of wave like equa- 
tions for three types of gravitational perturbations, 
according to the symmetry of the rotation group: scalar, 
vector and tensor. In four dimensions, scalar type of 
gravitational perturbations is called polar, vector type is 
called axial. Tensor type of gravitational perturbations is 
usually pure gauge in D = 4 black hole space-times. For 
D = 4 Reissner-Nordstrom- (anti)-de Sitter black holes, 
axial and polar types of perturbations are isospectral, so 
that for stability analysis, it is enough to analyze only 
one type of perturbations. For D > A, the isospectral- 
ity is broken, so that one needs to check all three kind 
of perturbations. The vector and tensor types of grav- 
itational perturbations were shown to be stable [3| for 
all D, with the help of the so-called S-deformation tech- 
nique. Therefore we shall consider here only the scalar 
type of gravitational perturbations. 

The equation of motion for gravitational perturbations 
of scalar type can be reduced to the wave-like equation 



i 



— +u;^-V±]^ir)=0, 



where the tortoise coordinate is defined as 

dr 



dr. 



fir)' 



V±ir) = fir) 



Mr^Hl 



(5) 



(6) 



(7) 



5m'' = 9^,u + ogf^v, 



(3) 



Here V+ and are potentials for the two kinds of 
scalar gravitational perturbations. The potential V_ re- 
duces to the pure gravitational perturbations when the 
black hole charge vanishes, while reduces to the per- 
turbations of the test Maxwell field in the black hole 
background in this limit. When the charge Q is non-zero, 
the gravitational and electromagnetic perturbations are 
coupled. 
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Note, that V+ is proven to be stable [3| with the help tential, which must be tested on stability, 
of the S-deformation. Thus we are left with the K_ po- Here we used the values 



i/-=A + ^&ll(l + A<5)X, 



(8) 



[/_ = [-A(f{d + 2){d+ 1)2(1 + \5fX^ + A8(f{d + l)id - 2)A(1 + XS)X 
-16{d - 2){d ~ 4)A2] Y - d^{M -2){d+ lfS{l + XSfX'^ 
-Ad'^id + 1)2(1 + XSf {{d + l)(3d - 2)XS -d'^}X^ 
+4{d + 1)(1 + XS) {X{d - 2){d - A){d + 1)(A + d^)5 

+Ad{2d^ - 3d + 4)A + d^{d~ 2){d - 4)(d + 1)} X"^ 
-16A {(d + 1)A (-4A + M^{d ~ 2)) 5 
+M{d - 4)A + U^{d + \)(d -2)]X 
+64A^ + 16d(d+2)A2. 

I 



(9) 



We shall imply that 

so that LOim > corresponds to a stable (decayed) mode, 
while LUirn < corresponds to an unstable (growing) 
mode. If the effective potential V{r) is positive definite 
everywhere outside the black hole event horizon, the dif- 
ferential operator 

h UJ^ 

drl 

is positive self-adjoint operator in the Hilbert space of 
the square integrable functions of r*, and, in that case 
any solution of the wave equation with compact support 
is bounded, what implies stability. An important fea- 
ture of the gravitational perturbations is that the effec- 
tive potential (EqE]), which governs the scalar type 
of the perturbations, has negative gap for the higher di- 
mensional black holes. Therefore the instability is not 
excluded for this case, and numerical analysis of pertur- 
bations is necessary. 
The values 



2X5 = 




l)2Af2 



1, 



A = (£-fd)(£- 1), £ = 2,3,4... 

are constants. 

For convenience we shall parameterize the black hole 
mass and charge by its event horizon and inner horizon 
r_ < r_|_ respectively. The value r_ = corresponds to 
the uncharged black hole. 



FIG. 2: Sample of instability for the Gauss-Bonnet black 
holes, D = 5, ^ = 2, a = 0.5. Instability corresponds to 
zero of Zi. 



III. NUMERICAL METHOD 

For analysis of stability we need to test a black hole 
response to external perturbations, which is dominated 
by the so-called quasinormal modes at late time. The 
quasi- normal boundary conditions correspond to the pure 
out-going waves at infinity and pure in-coming waves at 
the event (or de Sitter) horizon for asymptotically flat 
or de Sitter black holes. For asymptotically AdS black 
holes, the Dirichlet boundary conditions are imposed at 
infinity. If growing modes exist, the considered system is 
unstable. Although usually, damped quasinormal modes 
have both real and imaginary parts, i.e. are oscillating, 
the growing modes 8] are non- oscillating^ that is pure 
imaginary. This makes our search of unstable modes 
much easier. 

Below we shall discuss the numerical method, which we 
used here for asymptotically anti-de Sitter space-times. 

Let us start from the analysis of singularities of the 
equation ©. At the event horizon V±{r) cx /(r) ^ 0. 
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Therefore 



The quasinormal boundary conditions at the event hori- 
zon imply 

*(r) = (r-r+)-"^/^"('^+)(Zo + 0(r-r+)). 

At the spatial infinity the two linear independent so- 
lutions of \E'(r) are 

*i(r) ~r-(^-4)/2, *2W-r(^"6)/2, ^ 5,(10) 
^'i(r)-r~^/^ ^'2{r) r-^/^\n{r), D ^ 5, 

The quasi-normal boundary conditions imply that 

*(r oo) cx <( D ^ 5; (11) 

^-(D-4)/2 ^>6. 



Let us consider the new function 

i'^/f'{r+) 



2/(0 = 



«'(r). 



(12) 



If ^'(r) satisfies the quasi- normal boundary conditions, 
y(r) is regular at the event horizon. Since the function 
y(r) satisfies the linear equation, we fix its scale as 

yir+) = 1- 

Then ^'(^4-) can be found from the equation ([5]), 



where 



Vn 



2/'(r+)2 " (r+ - r_)/'(r+) + /'(r+) - 2iu 



Vq — lim 



r-^r\ f{r) 64r^iJ|(r+)' 

Imposing the above discussed boundary conditions at 
the event horizon, we solve the equation ^ numerically 
for each uj using the NDSolve built-in function in Math- 
ematica for r > rf, where 3> r+. 

In the general case the behavior of ^'(r) at infinity is a 
superposition of the two solutions (jlOp for in-going and 
out going waves $i(r), ^o{r), 



^/(r) -Z,$,(r) + Z<,$„(r), 



(13) 



where 5'o('') satisfies the quasi- normal boundary condi- 
tion (|lip . If w is the quasi-normal frequency, the corre- 
sponding solution must satisfy the boundary conditions 
(fTTj) at the spatial infinity and, thereby, Zi — 0. 

Thus, our numerical procedure is the following. We 
integrate the equation ^ numerically imposing quasi- 
normal boundary condition at the event horizon. At large 
distance we compare the obtained function ^'(r/) with 



(1131) and find, thereby, the coefficients Zi and Zq for any 
given value of w. The quasi- normal modes correspond to 
the roots of the equation 



Z^{lu) = 0. 



(14) 



In order to find Zi and Zo one has to find analytically 
expansions of $i(r) and ^o{r) at large distance. The 
expansion 



D = A] f d"^ 



r — —, D>5 
contains only inverse powers of r, while the expansion 



(15) 



1, D ^4 

r-^/'^\n(r), D ^5 
r^, D>6 




(16) 



contains also subdominant terms of the form of order 
ln(r) 

. Since the series in (jl6p are convergent we have 

r 

used only the first term of the expansion, which docs not 
contain logarithm. The expansion of p5p was done up 
to the order ~ r^^. 

The analytical expansion allows to find $i(r) and 
$o(r) within the desired precision for r ^ r^. li "^{vf) 
were known exactly, one would had found the coefficients 
Zi and Zo from the system of the linear equations 



*'(r/) 



*/(^/)- 



(17) 
(18) 



In practice, being the result of the numerical integration, 
the values '^(rf) and ^''(r^) contain a numerical error, 
which causes low precision of the coefficients Zi and Zo , 
found in this way. In order to minimize the numerical 
error, we find numerically the values of ^' at some large 
number of points near r — rf. Then we fit the obtained 
numerical values of ^ by the function f (r) . From 
the fit data we find Zi and Zo by solving the least squares 
problem at those points. 

Since unstable modes are purely imaginary, one can 
restrict the searching area for u by the imaginary axis. 
In this case the problem simplifies because the eigenfre- 
quencies and the coefficients Zi, Zq are real. It turns 
out, that the coefficient Zi changes its sign when cross- 
ing the solution. This can be used as an indicator of the 
existence of an unstable mode in the spectrum. 

In order to be sure that the above method indeed can 
find an instability we tested it for the two cases when the 
instability is determined both analytically and numer- 
ically by alternative method. Namely, we checked the 
instability of the black strings 01 1 01 1 and also found the 
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FIG. 3: Zi as a function of r_ and Im{uj) for D = 5, r+ = IR 




FIG. 4: as a function of r_ and Im{Lu) for D = 5, r+ = lOi? 



unstable modes of the Gauss-Bonnet black hole. Their 
values are in agreement with those, obtain within time- 
domain integration. Thus, as an example, on Fig. [2] one 
can see that the unstable mode is w = 0.18i, what per- 
fectly agrees with the value of unstable mode, obtained 
by the time-domain integration method in [loj |. 

Another restriction upon the possible values of unsta- 
ble modes comes from the depth of the negative potential 
gap {V — uo^ > Q guarantees stability), 

Im{uj) < \/-Vjmn, 

where Vmin is the minimal value of the effective potential 
at r-|_ < r < oo. 

For complex quasi-normal frequencies Zi is complex. 
For this case the solution of (fT4|) can be found by min- 
imizing Unfortunately, due to oscillation of the 
solution in the asymptotically flat and asymptotically de 
Sitter backgrounds, we were unable to fit the solution at 
very large distances. However, for asymptotically AdS 
background the solution does not oscillate at large dis- 
tance, and the described approach can be used to find 
quasi-normal modes of stable solutions. 




FIG. 5: Zi as a function of r_ and Im{u>) for D — 6, r+ = R. 




FIG. 6: Zi as a function of r_ and Im{uj) for D — 6, r+ = 
lOR. 

IV. STABILITY ANALYSIS 

For testing the stability, we have to perform the two 
tasks: to check that the there is no unstable modes by 
the method described in the previous section for the full 
range of the black hole parameters and, as a confirmation, 
to find the fundamental (damped, when the system is 
stable) quasinormal modes. 

We shall distinguish here three regimes: large black 
holes r-|_ ^ R, intermediate black holes r+ ^ R and 




FIG. 7: Zi as a function of r- and Im{u>) for D — 7, r+ — R. 
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FIG. 8: Zi as a function of r_ and Im{uj) for D = 7, r+ — 
WR. 




FIG. 9: Zi as a function of r_ and Im{uo) for D = 8, r+ = li?. 



small black holes <C i?, where i? is the anti-de Sitter 
radius. From the figures ([3|)- (fT6| . we can see that for 
the D — 5, 6 ... 11 black holes, Zi does not equal zero 
for any values of uj limited by Im{uj) < \J —Vmin,- We 
have shown this here mainly for the two values of r+i 
r_|_ ~ IR and r^. = 6R. These are representative cases of 
large and intermediate AdS black holes. For small black 
holes, an example of Zi behavior can be seen in fig. 1 171 for 
D — 5. There one can see that the smaller the size of the 
black hole, the larger Zi, what guarantees no-instability 
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FIG. 10: Zi as a function of r_ and Imiui) for D = 8, r+ = 

im. 




I I ( 1 . 11: Zi as a function of r_ and Ira(u) for D = 9, r+ = 
IR. 




FIG. 12: Zi as a function of r_ and Im{Lu) for D — 9, r+ = 
WR. 

at sufficiently small black hole size. Looking carefully 
at all range of parameters of r^, r_ and £, we have not 
found any zeros of Zi. Therefore we conclude that D = 
5, 6... 11 Reissner-Nordstrdm-anti-de Sitter black holes 
are stable for any values of the black hole parameters. 
Now, we shall check this by the search of the fundamental 
quasinormal modes, which, as it will be shown soon, all 
are damped. 

The spectrum of frequencies of neutral asymptotically 
AdS black holes is qualitatively different from asymptot- 
ically flat or de Sitter cases: the main striking feature of 
the spectrum is that almost all modes are proportional to 




FIG. 13: Zi as a function of r_ and Im{Lu) for D = 10, 
r+ = lO-R; logarithmic plot 
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FIG. 14: Zi as a function of r_ and Ira(u) for D 
r+ = IR; logarithmic plot 



= 10, 




FIG. 17: Zi for small Schwarzschild-AdS black holes D = 5; 
blue (top) for r-|- = R, red for r-f- = -R/2, yellow for r+ = 
green (bottom) for r+ = R/8. 
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FIG. 18: ReuJn=o for small D = 5 Schwarzschild-AdS black 
holes approaches the limit a; = 4 as r+ — » 



FIG. 15: Zi as a function of r_ and Ira{u) for D = 11, 
r+ = li? 



the radius of the black holes for large black holes r+ ^ R. 
The exception is the fundamental mode of the scalar type 
of gravitational perturbations of SAdS black holes: its 
real part approaches constant as r+ goes to infinity, while 
the imaginary part is inverse proportional to the radius 
of the black hole. Let us note that this property keeps 
also when AdS - black holes are charged. 

Quasinormal modes of a particular case of large D — h 
SAdS were considered in [201 ■ One can see in the table 




FIG. 16: Zi as a function of r_ and Im{Lu) for D = 11, 

r+ = im 



I that we accurately reproduce their results. Indeed, the 
Q = 0, r+ = 6 mode in our table I coincide with the 
fundamental mode of [l^ in proper units (see table III). 

Another property of asymptotically AdS black holes is 
that their quasinormal modes approach the real normal 
modes of the pure anti-de Sitter black holes, when the 
radius of the black hole goes to zero [l^ . This was shown 
for the test scalar fie ld p erturbations around D = 5 and 6 
SAdS black holes in [l^ . Here we have shown that for the 
scalar type of gravitational perturbations, quasinormal 
modes also reach their D > 6 pure anti-de Sitter values 

ojnR ^ 2n + D + £ - 3, D>6 pure AdS (19) 

The above formula for AdS space-time normal modes a;„ 
[2]| is valid only for D > 5, while for D = 5 the pure 
AdS spectrum is continuous, i.e. all modes are normal 
modes oi D = 5 AdS space-time. This happens because 
of the peculiar behavior of the effective potential at spa- 
tial infinity (which is zero in the tortoise coordinate): the 
effective potential has infinite negative pitch near r* = 
(see Fig. so that * ^ for all w [2l|. Notice that in 
spite of the infinite negative pitch, the area covered by 
the pitch is finite, so that the instability is not guaranteed 
a priori even for this case. The natural question arises: 
if the limit of small black holes has the same meaning 
for _D = 5 as it has for D > 6? On fig. [T51 one can see, 
that in the limit of small black holes, the D — 5 QNMs 
approach the same limit as D > 6 modes do. 

It is not remarkable that for Reissner-Nordstrom-AdS 
black holes, QNMs approach the same pure AdS values 
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TABLE I: Fundamental (n — 0) quasinormal modes of D = 5 
Reissner-Nordstrom-AdS black holes. 



r-/r+ 


D = 5, 


r+ = l 




D^5, 


r+ = 6 





2.20477 - 


0.58137 


i 


1.65171 - 


0.13759 i 


0.1 


2.19646 - 


0.57776 


i 


1.65143 - 


0.13623 i 


0.2 


2.17255 - 


0.56754 


i 


1.65058 - 


0.13214 i 


0.3 


2.13571 - 


0.55236 


i 


1.64914 - 


0.12544 i 


0.4 


2.08934 - 


0.53445 


i 


1.64713 - 


0.11639 i 


0.5 


2.03698 - 


0.51637 


i 


1.64455 - 


0.10549 i 


0.6 


1.98248 - 


0.50074 


i 


1.64150 - 


0.09350 i 


0.7 


1.93093 - 


0.48956 


i 


1.63800 - 


0.08131 i 


0.8 


1.88885 - 


0.48091 


i 


1.63404 - 


0.06989 i 


0.9 


1.85619 - 


0.46802 


i 


1.62944 - 


0.06059 i 


0.99 


1.89919- 


0.24619 


i 


1.62648 - 


0.05644 i 



TABLE II: Fundamental (n — 0) quasinormal modes of D = 
6, 7 Reissner-Nordstrom-AdS black holes. 



r-/r+ 


D = 6, r-+ = 6 


D^7, 


r+ = 6 





1.59784 


- 0.14933 i 


1.56442 - 


0.15517 i 


0.1 


1.59781 


- 0.14918 i 


1.56442 - 


0.15515 i 


0.2 


1.59761 


- 0.14814 i 


1.56438 - 


0.15492 i 


0.3 


1.59708 


- 0.14533 i 


1.56420 - 


0.15392 i 


0.4 


1.59603 


- 0.13994 i 


1.56372 - 


0.15123 i 


0.5 


1.59430 


- 0.13138 i 


1.56270 - 


0.14567 i 


0.6 


1.59175 


- 0.11954 i 


1.56086 - 


0.13607 i 


0.7 


1.58832 


- 0.10503 i 


1.55790 - 


0.12180 i 


0.8 


1.58395 


- 0.08923 i 


1.55359 - 


0.10364 i 


0.9 


1.57829 


- 0.07427 i 


1.54753 - 


0.08409 i 


0.99 


1.57267 


- 0.06664 i 


1.53968 - 


0.07200 i 



in the limit 7'+ (|19p . because one cannot assume 
r_|_ = 0, without taking Q 0. Thus, the fundamen- 
tal quasinormal modes of the scalar type of gravitational 
perturbations of D-dimensional RNAdS black holes obey 

u}nR^2n + D + £-3, r+->0, D>5 (20) 

Let us note that once we proved here the stability of the 
D-dimensional Schwarzschild-AdS black holes, the sta- 
bility of the Reissner-Nordstrom-AdS black holes can be 
intuitively understood from the behavior of the effective 
potentials (figures for scalar type in Q, 0]) at least for 
D > 6: the presence of the charge Q increases slightly 
the negative depth of the potential gap. Apparently the 
negative gap is not deep enough to allow bound states 
with negative "energy". The behavior of the effective 
potential for I? = 5 is quite different (fT7|) . yet, as we 
have shown, this does not lead to instability as well. 

V. DISCUSSIONS 

In this paper, by the numerical search of quasinor- 
mal modes, we have shown that Reissner-Nordstrom- 



anti-de Sitter black holes are gravitationally stable in 
Z? = 5, 6 ... 11 space-time dimensions. Before, the stabil- 
ity of asymptotically anti-de Sitter black holes was estab- 
lished only for D = 4 Reissner-Nordstrom-anti-de-Sitter 
black holes analytically 2]. Stability for D = 5 — 11 
found here and for = 4 found by Ishibashi and Ko- 
dama [2] does not contradict to the observed instability 
for D = 4 5 RNAdS black holes by Gubser and Mitra 
in [l^, [iJI, because the latter instability is induced by 
a tachyonic field coupled to the system in the Af = 8 
gauged supergravity. Thus, although metrics for the 
black hole in both cases are the same RNAdS metric, 
they are exact solution of different field equations, and 
the dynamic of perturbed equations is certainly different. 

The observed here stability of RNAdS black holes is 
interesting also, because we know that small AdS black 
holes {within the ordinary Einstein- Maxwell theory, con- 
sidered here) are thermodynamically unstable and may 
exert the Hawking-Page transition. It would be natu- 
ral to expect that this thermodynamic transition will be 
accompanied by a gravitational instability. Yet, as we 
have shown here, this does not take place, so that if the 
correlation between thermodynamic and gravitational in- 
stabilities exists, it is more subtle, than one could naively 
expect for complex gravitational systems. 

An important question, which was beyond the scope of 
our work, is the stability of extremally charged RNAdS 
black holes. Our closest aim is to give a detailed data on 
quasinormal modes of other types of gravitational pertur- 
bations (vector and tensor), and to find higher overtones 
of the spectrum j22| . We believe it would be interesting to 
investigate stability of charged asymptotically AdS black 
holes in the Gauss-Bonnet theory, where already there is 
instability, stipulated by Gauss-Bonnet terms, at higher 
multipoles. 
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